Frustration in c?-wave Superconducting Circuits: 7r-Ring Behaviour 
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A closed superconducting circuit containing an odd number of n -junctions, a 7r-ring, has a finite 
current in the ground state. We explicitly construct such rings for d-wave superconductors and 
demonstrate the existence of spontaneous currents by direct self-consistent solutions to the Bogoli- 
ubov de Gennes equations. We show that the current has a topological origin due to the frustration of 
the d-wave order parameter, which is only partially explained by the Sigrist-Rice tunneling formula. 
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It has very recently become clear that the flux quanta 
in a superconducting ring are in general of the form 
$ = (n + 7) < £ > o, where $ = h/2e, rather than simply 
$ = n$ as stated in the textbooks. The constant 7 
is analogous to the Maslov index in Bohr-Sommerfeld 
quantization for single particle dynamics. The predic- 
tion pL^I and subsequent experimental observation ||-||] 
of half-integer flux quantization (7=5) nas been one 
of the most dramatic developments in the field of super- 
conductivity in recent years. This discovery is profound 
for several reasons. The experiments provided clear and 
direct confirmation of the existence of d-wave supercon- 
ductivity in these materials. Previously, Wollman et al. 
had observed a closely related effect: split Fraunhofcr 
interference peaks in corner SQUIDs. However, the tri- 
crystal ring experiments ^-^| were the first to measure 
non-integral flux quantization directly. But the implica- 
tions of non-integral flux quantization go beyond simply 
confirming the pairing state symmetry, and have not been 
fully explored. The theoretical possibility of rings with 
flux quanta other than 7 = 1/2 or 7 = has been dis- 
cussed recently, f§ and this has not yet been observed. 
Rings with half-integer flux quantization also have some 
unique properties, including an exactly two-fold degen- 
erate ground state. This opens up new possibilities for 
novel superconducting devices, which may find applica- 
tions in quantum computing or elsewhere. The 
topology of half- integer quantization has been discussed 
recently by Volovik fll2|| . 

The original prediction of the half-integer flux quan- 
tization in d-wave superconducting rings was by Sigrist 
and Rice Q , following earlier work for the p-wave case . 
They considered a ring containing a single grain bound- 
ary junction, and assumed that the junction current was 
of the form 



I = I cos(2#i) cos(2#2) sin (p) 



(1) 



for d x 2_ y 2 pairing, where ip is the order parameter phase 
difference across the grain boundary junction, and 9\ and 
62 are the misalignment angles between the grain bound- 
ary normal and the crystal axes on either side of the 
junction. Provided that these angles are such that the 
junction is a 7r-junction (i.e. the pre-factor of sin(y>) is 



negative in Eq. [I]) and the ring inductance is large enough 
then the ring will exhibit a spontaneous current in the 
ground state. 

However, Eq. [I] probably does not apply. The grain 
boundary junctions have a complicated micro-structure 
and measured critical currents are approximately expo- 
nential in the misalignment angles |13[ . Also the current- 
phase relationship of the junctions is not simply given 
by sin (ip) juj. Furthermore the concept of a single 7T- 
j unction is not well defined since gauge transformations 
make ip ill defined for a single junction 0. Gauge invari- 
ant definitions can only be made for closed rings, and so 
one should talk about 7r-rings rather than 7r-j unctions. 

In this letter we show that the existence of 7r-rings 
stems directly from the frustration of the d-wave order 
parameter in certain ring topologies. The frustration can 
only be resolved by a spontaneous ground state super- 
current. Our calculations are the first to show this phe- 
nomenon at a microscopic quantum mechanical level. We 
show that for a frustrated d-wave ring there is a meta- 
stable energy maximum with zero current which has a 
zero energy state (ZES) at the Fermi level. In the ground 
state the ZES is removed, but at the cost of a sponta- 
neous circulating current. In contrast non- frustrated d- 
wave rings have no ground state current. Our numerical 
results for grain boundary junctions also differ consid- 
erably from the Sigrist-Rice model Eq. [lj; the current is 
not proportional to sin^j, and the angle dependence is 
not given by cos (Wi) cos (2#2). 

First let us review the symmetry principles underlying 
the Josephson characteristics of grain boundary junctions 
in d x 2_ y 2 superconductors. Assuming time reversal sym- 



metry, the current can be expanded |jl5fl 

I(6i,6 2 ,ip) = 7 ™ Wuh) sin (mp) . 



(2) 



For tetragonal lattice symmetry the functions I n (#1,6*2) 
are obviously periodic in 61 and 62 with period tt. Ro- 
tation of one crystal by 7r/2 changes the sign of the 
d-wave order parameter, and so we have a symmetry 
6*i — > 6*i +7r/2, ip — > iy9 + 7r. The implications of this sym- 
metry are different for current harmonics with n even or 
odd, 
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FIG. 1. a) The large angle grain boundary ring. The or- 
der parameter configuration shown is for 0° phase difference 
across the grain boundary, b) The [110] ring. The orientation 
of the order parameter lobes is for a phase difference of 0° 
across the [110] boundary. 

/»(0i + !,02) = -/n(0i,fc) nodd (3) 

In{0x + \,e 2 ) = +In{ex,6 2 ) neven. (4) 

The odd n harmonics have sign changes at certain angles, 
but the even n harmonics need not. The existence of tt- 
junctions is therefore not guaranteed by symmetry alone. 
Our earlier calculations |16| and experiments of IPichcv 
et al. |14j] show that the n > 1 harmonics cannot be 
ignored. These higher harmonics give a 4e//i Josephson 
coupling at angles where the 2e/h coupling vanishes [p"7| . 
In the weak coupling Josephson regime where only the 
n = 1 current component is significant the critical current 
I c = I\ (6i, 62) and this can be expanded as a Fourier 
series fl8| 



C cos 26»i cos 26> 2 + S sin 2B X sin 29 2 



(5) 



The Sigrist-Rice form, Eq. ED, is merely the first expan- 
sion term. If the other terms are significant then one 
cannot predict whether a given junction orientation has 
7r-junction or normal junction behaviour merely from the 
sign of cos 2^i cos 26*2 • Also, in contrast to Eq. |J the cur- 
rent does not vanish for the case Q\ — 82 = tt/4, corre- 
sponding to a twin boundary in a slightly orthorhombic 
crystal. 

In our calculations we directly compute the current in 
rf-wave superconducting circuits by solving the Hartree- 
Fock Gorkov or Bogoliubov-de Gennes (BdG) equations 
p9[ . We use a tight binding lattice with nearest neigh- 
bor hopping, t, and retarded attractive interaction LZy = 
— 3.5t with a frequency cutoff E c = 3.0t. We used a 
temperature of T = O.Oli and a chemical potential of 
/i = 0. The calculation is self-consistent in the non- 
local order-parameter, Ay = Uij{c^Cj\) , and in the hop- 



ping, tij = t + ^UijUij, where 



E 



,). Cur- 



rents are computed for each nearest neighbor bond [ p.6| , 
Iij, and self-consistency ensures that current conserva- 
tion is obeyed at each site. 
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FIG. 2. Evolution of order parameter phase within the 
ring. The arrow direction indicates increasing self-consistent 
iteration number so that by the time self-consistency is 
achieved we almost have a perfectly flat line (no phase gradi- 
ent) around the sample. 

To address the question of the origin of 7r-ring be- 
haviour at a microscopic quantum mechanical level we 
have directly calculated the current in rings of e?-wave 
superconductor. An example is shown in Fig. ^. This 
model grain boundary junction corresponds to misorien- 
tation angles of 61 = 2 = tan -1 1/2 = 26.6 deg. In ear- 
lier work we have studied the critical current and current- 
phase characteristics of such junctions in bulk supercon- 
ductors |l6)| . Here we have constructed superconducting 
rings containing these junctions. This is done by sim- 
ply connecting the edges labeled A and B in Fig. [J. 
The top and bottom edges shown in the figure cannot 
be connected, and so they are left with open boundary 
conditions. This provides a tractable microscopic model 
corresponding to the tri-crystal ring experiments of Tsuei 
and Kirtley |§-f§ . Of course the model has only one grain 
boundary, not three, but that is sufficient to demonstrate 
the principle. The other two junctions in the tri-crystal 
can be viewed simply as a means of connecting the two 
faces A and B. If the two additional junctions were iden- 
tical in symmetry then their own characteristics would 
be irrelevant |Q and the ring characteristics depend only 
on the single central boundary. In this ring geometry we 
have a system containing typically 200 or more inequiva- 
lent lattice sites and self-consistency in riij and A^- must 
be achieved independently on every bond. 

Fig. ^| shows the evolution of the phase of the ci-wave 
order parameter |l9| , |I(| in the ring system of Fig. |l| as 
the system approaches self-consistency. The phases plot- 
ted correspond to the sites indicated in black in Fig. [I]. 
Similar results are obtained for the other lines around the 
ring. The starting point for the calculation was (chosen 
arbitrarily) a step- function phase difference of tp = —120° 
applied across the grain boundary, with an opposite step 
at the A-B interface. Fig. || shows that as the itera- 
tive self-consistent procedure is carried out there is an 
'averaging-out' across the boundary, i.e the phase on the 
right hand side drops and that on the right hand side 



2 



FIG. 3. The [110] ring. The orientation of the order pa- 
rameter lobes is for a phase difference of 0° across the [110] 
boundary. 

increases until it reaches their mean value of —60° (a 
straight line) at self-consistency. Thus, for this system 
the ground state has no ground-state circulating current, 
and is therefore not a 7r-ring. 

The conclusion that this first geometry is not a Tr-ring 
is not surprising. In our earlier work J16| we calculated 
the supercurrent-phase characteristic, I(6i, 02, f) for a 
grain boundary junction between two bulk superconduc- 
tors for this geometry. Choosing a gauge so that zero 
phase difference corresponds to the d-wave orientation 
indicated in Fig. [l] leads to a positive current with a posi- 
tive phase difference, tp, as expected for a normal junction 
Jl6| . The sign of the d-wave order parameter can be cho- 
sen consistently around the circuit as indicated, implying 
the ring geometry of Fig. |l| has no intrinsic frustration 

Now contrast this with the ring geometry shown in 
Fiq. |^, where again edges A and B are joined. In this 
geometry the topology clearly indicates that the d-wave 
order parameter is frustrated |L5| . Suppose we define the 
gauge so that zero phase difference across the junction 
corresponds to that indicated in Fiq. |^. Then one can 
see that one must introduce a sign change at the point 
where the A and B edges are joined to make a close 
ring. Alternatively in a different gauge convention one 
could impose no sign change at the A-B interface, but 
then necessarily there must be a sign change elsewhere 
in the loop. This intrinsic frustration is a new feature of 
the system Fiq. [| and is quite different from Fig. |l| where 
there was no similar frustration. Topologically Fiq. ^ is a 
globally non-orientable surface, equivalent to the Mobius 
strip (one cannot continuously assign a unique orienta- 
tion (gauge) at each point), while Fiq. [l] is globally ori- 
cntable and topologically equivalent to the surface of an 
ordinary cylinder. See Ref. fl^| for further implications 
of this topology. 

The dramatic effect of this topological difference is ap- 
parent in the evolution to self-consistency, as seen in 
Fig. [|. Starting the evolution to self-consistency with 
Aij = for all sites on the left hand side of the bound- 
ary and a bulk d-wave state with an arbitrary phase of 
30° on the right, we obtain the order parameter phase 
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FIG. 4. Evolution of the phase difference across the [110] 
interface with increasing self-consistent iteration number. 
The arrows point in the direction of increasing iteration num- 
ber. 

shown in Fig. ^|. The phase becomes increasingly lin- 
ear, corresponding to a uniform and non-zero supercur- 
rent throughout the ring. The system has therefore a 
spontaneous ground state current and is therefore a Tr- 
ring. Note that the discontinuity in phase at the join 
of the A and B edges is purely due to the sign con- 
vention chosen in Fig. [| and does not correspond to 
any physical discontinuity in the Ay . Also the slight 
non-linearity in the phase gradient at x = is consis- 
tent with the increased cross-section width of the strip 
at this point. Fig. |]is a direct microscopic example of a 
Tr-ring due to topological frustration in d-wave supercon- 
ductors. In order to understand the microscopic origin of 
this spontaneous current it is helpful to first consider a 
single 9i = 62 = 45 deg junction between two bulk super- 
conductors. This would be a [110] twin boundary for two 
slightly orthorhombic superconductors. The geometry is 
identical to Fig. || except that the edges A and B are not 
joined, but rather connect to infinite bulk superconduc- 
tor on either side of the junction. Now there is transla- 
tional symmetry parallel to the interface, and every site 
is symmetrically equivalent to one of those indicated in 
black on Fig. ^. In this case we can use the numerical 
methods of Ref. Jl6| to determine the characteristics of 
this junction. We choose to define zero phase difference 
ip as indicated in Fig. [|. In this case the non-local order 
parameter A y - clearly has a perfect bulk d-wave shape at 
every site, including the sites at the center of the bound- 
ary. Calculating the local density of states for this does 
indeed yield an ideal bulk d-wave density of states. 

Suppose now that we apply a phase difference tp across 
the interface. Rotating the order parameter on the right 
hand side of the interface by 90° corresponds to a phase 
difference of <p — 180°. In this case the local density of 
states is shown in Fig. ^|(a). There is a strong zero energy 
state (ZES), or resonance, at the center of the gap. This 
ZES is essentially of the same origin as the one described 
by Belzig, Bruder and Sigrist pi] . 
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FIG. 5. a) LDOS for the [110] interface in bulk with a phase 
difference of 180°, b) Current versus phase for same boundary, 
c) 0° phase difference across the [110] interface when joined 
into a ring, d) [A<j| for same system as c): full line is for 30°, 
broken line for 0°. 

However in our case we do not observe splitting of the 
peak, and there is no time reversal symmetry breaking 
(TRSB). The final self-consistent answer yields a d-wave 
order parameter but with a small extended-s component 
at this phase difference. The pairing state at the junc- 
tion is locally s + d and not s + id. We have also calcu- 
lated the supercurrent-phase characteristic of this [110] 
twin boundary junction. We apply a phase difference tp 
between the two bulk superconductors and calculating 
self-consistently the order parameter Ay , charge ny and 
current Zy on each bond in the region of the junction. 
Fig.^(b) shows that the current versus phase profile is an 
almost linear saw-tooth function with sharp discontinu- 
ities at ±7r. It is qualitatively the same as for the other 
boundaries previously reported 16 1. The shape is indica- 
tive of a strong coupling junction [p2[ , and we attribute 
the sharp discontinuity to the existence of the ZES at 
ip = 7r. Obviously the fact that the curve is not simply 
a sine wave indicates that the higher harmonics are rele- 
vant in Eq.^|, consistent with the experiments of Il'ichev 
et al. juj. Also note that in this junction geometry the 
Sigrist- Rice formula Eq. |l] predicts no supercurrent flow 
at all, but the non-zero current is due to the sin 26\ sin 2#2 
or higher order terms in the expansion Eq. || 

The state ip = n in Fig. ||(a) with the ZES is the en- 
ergy maximum for the junction, as can be seen by inte- 
grating the current, J I dip, of Fig ||(b). It corresponds 
to a phase slip of tt in the order parameter across the 
junction. Such states are also possible in the ring geom- 
etry. Fig. H(c) shows a self-consistent order parameter 
phase for the same ring geometry of Fig. [|. This was 
obtained by starting the self-consistent calculation with 
a bulk e?-wave order parameter and with the phase ar- 
rangement shown in Fig. |^ Thus, the calculation begins 
with a phase-slip at the A-B boundary. At this point the 
magnitude of the <i-wave order parameter |A,j| is heavily 



suppressed, as shown in Fig. ||(d). This is because locally 
there are sites where the Ay to the four neighboring sites 
are (+, — , — , +) rather than (+, — , +, — ) counting clock- 
wise. We interpret this phase slip solution for the ring 
as the meta-stable energy maximum state which sepa- 
rates the two equivalent energy minima corresponding to 
positive or negative circulating current. 

In summary we identify ring geometries in c?-wave su- 
perconductors as either frustrated or not. Frustrated 
rings have a circulating current in the ground state for 
but no ZES. The only available states without current 
flow have a phase slip, but with the onset of a non-split 
ZES (i.e. no time reversal symmetry breaking). We also 
calculated the current versus phase profile and demon- 
strated that the Sigrist Rice formula is not a valid ap- 
proximation. Finally, we suggest that experiments de- 
signed to measure the differential conductance in these 
rings are an important next step in understanding the 
role of the ZES. 
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